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Abstract — We present a new algorithm, A*Pnine, to list (in order of in- 
creasing length) the first K Multiple-Constramed-Shortest-ftth (KMCSP) 
between a given pair of nodes in a digraph in which each arc is associated 
with multipie Quality-of-Service (QoS) metrics. The algorithm constructs 
paths starting at the source and going towards the destination. But, at each 
iteration, the algorithm gets rid of all paths that are guaranteed to violate 
the constraints, thereby keeping only those partial paths that have the po- 
tential to be turned into feasible paths, from which the optimal paths are 
drawn. The choice of which path to be extended first and which path can 
be pruned depend upon a projected path cost function, which is obtained 
by adding the cost already incurred to get to an intermediate node to an ad- 
missible cost to go the remaining distance to the destination. The Dykstra's 
shortest path algorithm is a good choice to give a good admissible cost. Ex- 
perimental results show that A*Prune is comparable to the current best 
known e-approximate algorithms for most of randomly generated graphs. 
BA^Prune, which combines the A^Prune with any known polynomial time 
e-approximate algorithms to give either optimal or e-approximate solutions 
to the KMCSP problem, is also presented. 

Kepvords — shortest paths, constraint based routing, QoS routing, multi- 
ple constrained path selection, Dgkstra algorithm, NP complete. 

I. INTRODUCTION 

QUALITY of Service (QoS) sensitive routing is of critical 
importance in achieving a data network with high speed 
and efficiency (high ratio of capacity/cost), as well as meeting 
QoS for each connection. While efficient utilization of network 
resources is important to the service provider, the user is in- 
terested in both the speed and the quality of service provided. 
The need for QoS routing can be justified for both reservation- 
based services (e.g., Intserv, ATM) as well as reservationless 
services (e.g., Diffserv). The goal of QoS routing is to find 
a path that satisfies multiple QoS constraints while achieving 
overall network resource efficiency [1]. A key challenge in QoS 
routing is the K-Multiple-Constrained-Shortest-Path (KMCSP) 
problem, which is to find K feasible Constrained-Shortest-Path 
(CSP) from a source node to a target node and subject to mul- 
tiple constraints (e.g., bandwidth, delay, jitter, administrative 
weight), and list them in order of increasing length. One special 
case of KMCSP problem is the Multiple-Constrained-Shortest- 
Path (MCSP), which is the KMCSP problem with K = 1, 

The MCSP problem is a basic optimization problem that is 
both theoretically interesting and have many practical applica- 
tions [2], and hence has received considerable attention in the 
literature [3] - [8]. 

Besides all the applications that are applicable to the MCSP 
problem, the KMCSP problem can also have many other prac- 
tical applications, such as in the network routing and network 
design problems. 
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In general, network routing consists of two basic tasks [5]: 
distributing the network information and searching for best fea- 
sible path with respect to given constraints. We focus on the 
second task and assume that the network information is avail- 
able to every node (e.g., via link-state routing). Each link in the 
network is associated with multiple QoS metrics, which can be 
either static or dynamic (based on whether varying with time), 
either additive or non-additive (based on the additivity along a 
path) [7]. For the additive parameters (e.g., delay, hops, jitter), 
the cost of an end-to-end path is given, exactly or approximately, 
by the sum of the individual link values along that path. In con- 
trast, the cost of a path with respect to a non-additive parameter, 
such as bandwidth, is determined by the value of that constraint 
at the bottleneck link. In the case of dynamic routing, the band- 
width is a dynamic parameter since its value may change each 
time a demand is routed or torn down, while the length of a 
link is a static parameter since it remains unchanged. In this pa- 
per we will mainly focus on additive and static parameters. The 
constraints associated with non-additive or dynamic parameters, 
can be dealt with first solving the KMCSP problem subject to all 
the additive and static constraints, and then select one path from 
these K CSPs such that all the other non-additive or dynamic 
parameters are satisfied. Most of the current approaches handle 
the non-additive or dynamic constraints by first pruning out all 
links that do not satisfy these constraints and then solving the 
MCSP problem in the residual network. The MCSP procedure 
need to be called each time a new demand need to be routed, 
and it may take time to find a feasible route if the MCSP pro- 
cedure is time consuming, while our KMCSP method can speed 
up the online routing time, since selecting a feasible path from 
the K precomputed candidate paths is generally much faster than 
solving a MCSP problem in a residual network. 

Generally, a network design problem is to find a least cost 
or a maximum revenue network, such that a given set of de- 
mands are routed through routes which meet some given QoS 
constraints. The KMCSP algorithms can be used to pre-compute 
a set of the candidate paths, then the network design problem 
can be formulated as a path based Linear Programming (LP) or 
Integer Programming (IP) problem, or can be solved by select- 
ing paths from the precomputed path Hst through some other 
heuristic methods. We have appHed the KMCSP algorithm in 
SPIDER, a web based optical network design tool [9]. 

In the case of single metric, the KMCSP problem becomes 
a problem of finding if -Shortest-Paths (KSP) from one source 
node to one target node. The KSP problem has many practi- 
cal applications [10] [11]. In multiple constrained metrics case, 
both the MCSP and the KMCSP problems are known to be NP- 
complete[3][4][12][14]. 

Most of the current approaches are concentrated on develop- 
ing efficient polynomial or pseudo-polynomial-time algorithms 
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to give feasible or approximate solutions to MCSP problem. Ex- 
isting algorithms can be grouped into three approaches: 

1. e-approximate solutions given by pseudo-polynomial-time 
algorithms in which the complexity depends on the actual val- 
ues of the link weights (e.g., maximum link weight) in addition 
to the size of the network [4] [14]. 

2. Polynomial-time heuristics that are fast in searching for a 
feasible path but cannot guarantee finding one in the network 
[2][3][5]. 

3. e-approximate algorithms that are polynomial [12][13]. 
The main contribution of this paper is developing an algo- 
rithm to give exact solutions to the KMCSP problem, and can 
also be applied to its special cases such as MCSP and KSP prob- 
lems. The algorithms presented in this paper can be applied to 
multiple constraints and can find any required number of CSPs. 
To give an exact solution, the running time may be exponential 
in the worst case. However, our test cases show that the actual 
running time of our algorithm is comparable to the existing ap- 
proximate algorithms for most of practical networks. 

The rest of the paper is organized as follows. The KM- 
CSP problem and its solutions are further analyzed in Section 
n. A*Prune, an algorithm for KMCSP problem, is presented 
in Section III. Section IV analyzes the performance of the 
A*Prune algorithm. Section V discusses some efficient meth- 
ods for computing the lower-bounds, which are used to give 
the look-ahead feature in A*Prune, thus variant of A*Prune 
corresponding different lower bounds for the KMCSP, such as 
A*Uniform, A*Dijkstra, are presented. In Section VI, we com- 
bine the A*Prune with any known e-approximate algorithm to 
give the Bounded A*Prune (BA*Prune), which can give ei- 
ther exact or e-approximate solutions to KMCSP problem in 
polynomial-time. In Section VII we provide some experimental 
results and the comparisons for the algorithms presented in this 
paper and some well-known algorithms. Conclusions are pre- 
sented in Section VIII. All the proofs to the lemmas are given in 
the Appendix. 

II. Problem Formulation 

The KMCSP problem can be defined as follows. 

Definition 1: KMCSP problem: Consider a network that is 
represented by a graph G = {V^ E), where V is the set of nodes 
and E is the set of links. Each link e E is associated 
with R non-negative and additive QoS values: Wr{i,j),r = 
1, 2, . . . , iZ. A length (sometime called cost) function wq is de- 
fined as follows: 

R 

yjoihj) = Yl^rWrihj) (1) 
r=l 

Given a source node s and a target node t, and jR constraints 
Cr (s, t), r = 1, 2, . . . , R The KMCSP problem is to find either 
the first K shortest length paths or all the paths (depending on 
which number is smaller) from a source node 5 to a target node 
t such that 

»Vr(p(s,t)) = E Wr{i,j)<Cr{s,t) 

(<.j)6p(»,t) (2) 

Vr€(l,2,...,/J). 



We use a bold character to represent an {R + l)-dimensional 
vector in this paper, such as X E represents X = 

{Xo, Xi , . . . , Xji), Generally, the constraints vector C(z, j) de- 
pends on the node i and node j, however, we sometimes use C 
to represent C(s,*). 

In the case of iiT = 1, the KMCSP, or simply the MCSP prob- 
lem can be formulated as an integer programming (IP) problem 
as follows. 

Definition 2: MCSP (KMCSP with K = 1) as an (edge 
based) IP problem : 

given : 

G(y,E),se v,t e v,r> i,k > i; 

wr{ij)>0, V(i,j)e£?, Vr€(l,2,...,ii); 
ar>0, Cr>0, Vr € (l,2,...,iJ); 
define : 

R 

'^oihj) = E 0,rWr{iJ); 

Wr^ E'wriiJ)x{iJ),^re{0,h...,R); ... 

minimize : Wo; 
subject to : 

a:(i,j)e{0,l}, V(i,j)€E; 
E y^r{iJ)xiiJ)<Cr, Vre(l,2,...,ii); 

f 1, if i=S 

E ^ihj)- E a:0'>0 = S -l> »/ » = <^ 
i:|(i,i)e£? j'AU.i)eE [ 0, else 

Let dip) be the destination node of the path p in this paper. 

To cope with the KMCSP problem, we give notations to some 

sets of paths as follows. 

P{V, V) = {x: \x is a path in G}; 
P(i^ V) = {x :\x e P{V, V) and i is its source node}; 
P{hj) = {x:\x ^ P{iy V) and j is its target node}; 
P(z,K,f(x),C(2,f)) = 

{x:\xGP{i,V),f{x)<C{i,t)}; (4) 
P(z,j,f(x),C(i,j)) = 

{x:\xe P{i. V, f (x), C(2, j)) and d{x) = j}; 
P{i,j,{{x),C{iJ)yK) = {x :\x e the first K shortest 

length paths of P(i, j, f (x), C(i, j))}; 

Here, Co = oo, i e V, j e V, and f (p) G E^-^^ repre- 
sents (R+1) given functions of the path p. The constraint Co is 
redundant and thus can be set to infinity. 

Using these symbols, the solution sets of different problems 
can be represented as follows. 

KMCSP: P(5,t,W(j>),C,Jf) 

MCSP: P{s,t,W{p),C,K)\K=i (5) 
KSP: P{s,t,W{p),C,K)\n=i,c^=oo 

Definition 3: Simple path: A simple path is a path without 
loops. 

Definition 4: Head path and tail path: Let node n be an 
intermediate node of the pathp(i, j). Node u divides pathp(z, j) 
into two paths, path p(i, u) and path p(u, j). then path p(i, it) is 
called a head path of path p{i, j), path p(w, j) is called a tail path 
of path p(ij j). We represent the path p{ijj) as the combination 
of a head path and a tail path in the following format: 

P{hj) =P{hy)p{uJ) (6) 
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Definitions: Additive Parameters: Let Wr{p(iJ)) be a 
parameter associated with the path j), and Wr(p{i,j)) = 
Wr{i, j) when p{i, j) is only one link path, then Wr is called an 
additive parameter if 

Wr{p{i,u)p{uJ)) = Wr(p{hu)) + WripiuJ)) (7) 

remains true for all i E Vyj e V,u eV. 

Definition 6: lower-bound distance: Let W(p(i,j)) be a 
vector associated with the pathp(z, j) and as defined in (2), then 
a vector D(i, j, C(i, j)) is called a lower-bound distance vector 
from node i to node j and associated with the constraint vector 
C{iJ) if 

Dr{i, j, C(i, j)) = min Wr(p{h j)) 

Vr 6(0,1,. 

Definition 7: Admissible distance: A vector A(i, j, C(i, j)) 
is called an admissible distance from node i to node j and asso- 
ciated with the constraint vector C(i, j) if 

0 < A(i, j,C(i,i)) < D(z, j,C(i, j)) (9) 

Here, D is a lower-bound distance and the operator < with two 
vector operands is defined as: 

SL<h^ ar <br, Vr G (0,l,...,i?) (10) 
Definitions: Projected distance: Given a path p(s,i) € 
P(s, V), a constraint vector C(5, and an admissible distance 
A(i, j, C(2, j)), a path function H(p(5,i)) is called a projected 
distance associated with the node pair (s, t) and the constraint 
C(5,t) if 

H(p(5,i)) = W(p(5,i))+A(t,t,C(s,t)-.W(p(s,i))) (11) 
Definition 9: Feasible path: A path p{s, t) is called a feasible 
path associated to the node pair (s,f) and the constraint vector 
C{s,t) ifp{s,t) £ P(s,f, W(p), C(s, t)) Within this meaning, 
P(s, t, W(p), C) can also be called the feasible path set 

Definition 10: admissible head path: A path p(s, i) is called 
an admissible head path associated to the node pair (s,t) and 
the constraint vector C(5,i) ifp(s,i) E P(s, y,H(p(s,i)), C). 
Here H(p(s,i)) is a projected distance. Within this meaning, 
P(i,F,H(p(s,i)),C) can also be called the admissible head 
path set. 

Based on the definitions given above, we make the following 
observations. 

Umma 1: P{s, t, H(p), C, K) is the solution set of KMCSP 
problem, i.e., 

P{s,t,H{p),C,K)=P{s,t,W(p),C,K) (12) 

Lemma 2: Vp 6 P(5,V), path p can always be expanded 
from the trivial path p(s, s). 

Lemma 3: A path expanded from an inadmissible head path 
must be an inadmissible head path, thus can never be a solution 
path, i.e., 

q ^ P{s, y, H(p), C) 92/ P(5, y, H(p), C) (13) 

\/qePiV,V) and \/y e P{V,V) 
These Lemmas tell us that we can get the solution set of KM- 
CSP by expanding the trivial path p(s, s) and all its extended 
feasible head paths step by step. This gives the basic ideas of 
the A*Prune Algorithm. 



III. A*Prune Algorithm 

A*-search[15][16][17][18], as well as uniform-search, 
breadth-first-search and depth-first-search are well known 
searching strategies in Artificial Intelligence [19]. We com- 
bine the A* -search with a proper pruning technique to get the 
A*Prune algorithm, which can be used to solve the KMCSP 
problem. Using the terminologies described in the previous sec- 
tions, the A*Prune algorithm can be simply described as: start- 
ing from expanding the path p{s, $), potentially, all the paths in 
P(s, V) can be reached; however, with a proper pruning against 
the given constraints C, only the paths in admissible head path 
set P(5, Vy H(p), C) remain as candidate paths for further ex- 
panding; furthermore, the candidate paths are ordered properly, 
such that the path with shortest projected length Ho(p) is se- 
lected and expanded first, then we can terminate our expansion 
procedure once we have found enough number of CSPs or there 
are no candidate paths left. This procedure will only expand a 
subset of the admissible head path set P(s, V, H(p), C). 

A pseudo-code of A*Prune algorithm is shown in Figure 1. 

The key processes in A*Pnine are explained and analyzed as 
follows. 

1. Pre-compute an admissible distance, A(i, ^, C(z, *)), Vi G 
V: Here C(i,t) represents some known constraints associated 
to the node pair (i,t). We can set C(i,t) = oo if nothing is 
known about C{i,t). Many existing algorithms, such as Di- 
jkstra*s shortest path algorithm can be used in finding a good 
admissible distance. We will look at the question of designing 
good admissible distances in Section V. 

2. Path expanding: Suppose we have an admissible head path 
list, AHP-heap, which is initialized to contain the trivial path 
p{sys). The path expanding process first selects and removes 
a path p from AH P .heap, then expands the selected path one 
step further to get all possible extended paths and inserts all the 
admissible head paths into the AHPJieap. 

3. CSP collecting: Each time we take a path from AHPJkeap, 
we first check if the path is a CSP. If it is, the path is saved as a 
solution and will not be expanded further. 

4. Candidate path list ordering: If the candidate path list 
AHPJieap is ordered in a way such that the path with shortest 
projected length Ho{p), and the largest length Wo{p) if breaking 
the tie, is put on the head of AH P -heap, then we can stop the 
expanding process once we have found enough number of CSPs 
without any loss of the optimality. Once the candidate list is 
ordered in this way, the A*Prune search will be best (with short- 
est Hoip)) first search, and depth ( largest Wo(p)) first search if 
several paths has the same value of Hoip)- The heap-sort[20] 
algorithm is used here for its efficiency. The candidate path 
list AHPJieap is heap-sorted whenever a path is added or re- 
moved. 

5. Inadmissible head path pruning: Once a new path is gen- 
erated in the path expanding process, a check is made against 
the given constraints, using some lookahead features. The 
newly generated path will be put into the candidate path 
list AHPJieap if it belongs to the feasible head path set 
P(s, V^,H(p), C). All the inadmissible head paths are pruned 
out and will not be expanded further. So, if the path ex- 
panding process is combined with the constraint pruning pro- 
cess, all the paths in AHPJieap are admissible head paths. 
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function A*prune{G, 5, t, w(jE7), C, K, R) 
inputs: 

G = (V; jE7), a graph with node set V and edge set E; 

{s, t): a node pair with source s and target t\ 

K: number of paths to be found; 

R: number of constraints, R= I 'lf applied to KSP; 

w(e): R metrics associated to each link eS E; 

C = C{syt): R constraints; 

1. Vi G V, 2 7^ 5, i t and Vr e (1, 2, . . . , R), compute: 

Dr{iy t); H length of Dijkstra path from i to ^ 
Dr{s^ %)] II length of Dijkstra path from s to i 

Ar(i,t,C(i,*)) ^ -Dr(i,t);//Addmissible distance 

2. Initialize: 

ik ^ 0 ; // number of CSPs found 

W(p(5,5))^0; 

H(p(5,s))^ A(5,t,C(s,t)); 

AHPJieap i- {p($, s))\ll admissible head path heap 

CSPJist ^ { } ; //set of found CSPs 

3. wiiile (k <K and AHPJieap is not empty) 



4. g(s, u) ^ the first path of AHP.heap; 

5. Remove the first path of AHPJieap; 

6. heapsort{AHP.heap) ; 

7. w 4- the end node of u); 

8. ifw = Uhen 

9. insert q{s, u) into CSPJist; 

10. k^k + 1; 

11. goto 3; 

12. end if // q(s, u) is saved as a CSP solution 

13. OutEdges -<-{all edges outgoing node u}; 

14. while (OutEdges not empty) 

15. e (u , t;) ^ a removed edge from OutEdges ; 

16. p(s, v) q{$, u)e{u, v); 

17. W(p(5, v)) ^ W(g'(5, u)) + w(e(u, v)); 

18. if simple path is required and v E ^(s, u) 

19. then goto 14; 

20. end if //non-simple head paths are pruned 

21. forr = l,2,...,iJ 

22. if Hr(p{s, v)) > Cr then goto 14; 

23. end if 

24. end for // inadmissible head paths are pruned 

25. Insert v) into AH P. heap; 

26. heapsort{AHP-heap) ; 

27. end wiiile // p{s, u) has been expanded 



28. end while 

29. return CSPJist; 

30. stop 

function Hr[p{s, v)) 

reinms Wt{p{s,v)) ^- Ar{v,tyC{v,t)); 

function heapsort{AHPJieap) returns a heap 

AHPJieap such that the path p(s,i) with minimum 

^o(p(5,i)), and maximum Wo(p(s,0) case of breaking 

the tie, is put on the head of the heap AHPJieap. 

Fig. 1. Algorithm of A*Pruiie 



0-7803-7016-3/01/$10.00 ©2001 IEEE 



i.e., AHPJieap C P(s, V,H(p),C). Furthermore, all the 
paths in P{s, V, H(p), C), will be eventually expanded once the 
AHPJieap is exhausted to empty. 

6. non-simple path pruning { required for simple path searching 
only): This step is required only if simple path is an additional 
requirement to the CSPs. Any non-simple path does not need to 
be put on the candidate path list, since any of its extended paths 
must also be non-simple. Suppose all the paths in the candidate 
path list are simple paths, then a newly generated path extended 
from any of the candidate paths is simple if its end node presents 
only once in the path. So we can prune out all the non-simple 
paths by checking the appearance of the end node in a newly 
generated path. 

7. Terminating condition: The A*Prune program will be termi- 
nated either it has found the required number of CSPs or there 
are no paths left in the candidate path list. The candidate path 
ordering steps let the candidate path with shortest projected path 
length, and the longest length path if breaking the tie, to be ex- 
panded first. So the ordering steps may result in finding the 
required number of CSPs as early as possible before exhausting 
all the paths in AHPJieap. The pruning steps try to keep the 
AHPJieap to contain candidate paths as less as possible. 

The A*Prune algorithm combines all these processes to se- 
lect, expand, prune the candidate path list step by step, until the 
required number of CSPs are found or there are no candidate 
paths left. 

IV. The Performance of A*Prune 

The performance of an algorithm is usually evaluated in terms 
of the following four criteria [19]: 

• Completeness: is the algorithm guaranteed to find a solution 
when there is one? 

• Optimality: does the algorithm find the highest-quality solu- 
tion when there are several different solutions? 

• Time complexity: how long does it take to find the solution? 

• Space complexity: how much memory does it need? 
Lemma 4: The completeness of A'^P rune: A*Prune is com- 
plete in finding simple KMCSP paths. However, A*Prune may 
not be complete in finding complex KMCSP paths. 

Lemma 5: The optimality of A *Prune: A*Prune is optimal in 
finding either simple or complex KMCSP paths. 

Let Q be the number of expanded paths, then, the space com- 
plexity of A*Prune is Q, and the time complexity of A*Prune 
is dQ{R + h + log Q), Where K is the number of paths to be 
found, h is the maximum hops of these K shortest paths, d is 
the degree of G and R is the number of the constrained metrics. 
However, A*Prune can still be exponential, since Q may be ex- 
ponential. It has been proved that A* can be sub-exponential 
growth if the error in the heuristic function grows no faster than 
the logarithm of the actual path cost[21][22]. Because of the 
pruning, our A*Prune is more efficient than A*. Thus, the above 
conclusion for A* also remains true for A*Prune. In mathemat- 
ical notation, the condition for sub-exponential growth is that 

|D(z, t, C(i, t)) - A(i, t, C(2, t))\< 0(log |D(i, t, C{i, t))\) 

Therefore, it is always better to use an admissible distance 
A(i, t, C(i, t)) with higher values, as long as it does not overes- 
timate the lower-bound distance D(t, t, C(i, t)). 

3 IEEEINFOCOM2001 



When applied to KSP problem, A*Prune becomes polyno- 
mial growth algorithm and the length of candidate path list is 
bounded by: 

Q<Khd<KN^ (14) 

Then the space complexity of A*Prune for KSP is Khd and the 
time complexity is 0{Kh(f(h + ii + log{Khd))), 

V. Methods for Computing Admissible Distance 

So far we have seen the two extreme ends of admissible dis- 
tances: A(i,t, C) = 0 and 'D{i,t,C{i,t)) given in ( 8). The 
former one does not need any computation but does little help 
to the look ahead search, while the second one is really helpful 
to the look ahead search but it also increases the precomputa- 
tion time. We need some admissible distances between the two 
extreme ends, such that they are close to the best admissible dis- 
tance and can be computed efficiently. 

A problem with less restrictions than the original problem is 
called a relaxed problem. It can be easily shown that a solu- 
tion of any relaxed problem of the MCSP problem, as shown 
in (3), always gives an admissible distance A{s,t, C(s,t)). In 
this section, we present some efficient methods for computing 
admissible distances by solving a relaxed problem of (3). 

A. Uniform Distance 

Obviously, A(i,t, C{i,t)) = 0, Vi G V, is an admissible dis- 
tance, which is called uniform distance. A*Prune with uniform 
distance is called A*Uniforni algorithm. 

A*Uniform, like A*Prune, expands all paths from source 
node. However, instead of using heuristic functions as look- 
ahead feature, A*Uniform checks the actual path criteria for 
pruning and the shortest length path is expanded first. 

A*Uniform seems too greedy and probably is not as efficient 
as A*Prune with greater admissible distance. However, in the 
case of finding CSPs from one source node to N target nodes 
and with uniform constraints (C(2,t) = C{s,t), Vi € V), 
A*Uniform maybe more efficient than running N rounds of 
A*Prune. 

B, Dijkstra Distance 

For each metric Wr, we can use Dijkstra 's algorithm to find 
the path with shortest Wr(p). We can find R shortest metrics 
associated with each r e (1, 2, . . . , jR). These R shortest met- 
rics make up a vector, which is called the Dijkstra distance and is 
represented by the symbol D(i, oo). By definition and using 
(8), we have: 

Drihhoo) = min Wr(p{hj)) (15) 

Vre(0,l,...,iZ) 

Lemma 6: The Dijkstra distance is admissible. 

Therefore, A*Prune can use Dijkstra distance as an admissi- 
ble distance and is called A*Dijkstra algorithm in that case. The 
A*Dijkstra algorithm remaines unchanged as shown in figure 1 
except the procedure of updating the admissible distance at lines 
?? and ?? can be omitted. 



In the case of simple path is also a requirement to KM- 
CSP problem, we can compute the Dijkstra distance in network 
G' = G - s (the network G with node s and all the edges linked 
to node s removed), since we really do not want any of the looka- 
head paths passing through node s again. Usually, the Dijkstra 
distance in network G' can give a better admissible distance than 
the Dijkstra distance in network G. 

C. Admissible Distance Given by Other Known Heuristic Algo- 
rithms 

There are many known algorithms which can give lower 
bounds to the MCSP problem [12]. Given an e-approximate 
solutions [14], a strict lower bound can also be retrieved by sub- 
tracting the error from the approximate solution. Lagrangian 
relaxation method [23] can also be used to generate a lower- 
bound to the IP problem listed in (3). Obviously, any of these 
algorithms can be used for computing the admissible distance in 
our A*Prune. 

VI. BA*Prune 

Despite all the nice properties of the A*Prune algorithm, such 
as optimality, completeness, ability to solve KMCSP, A*Prune 
is not polynomial time. Both the time and space complexities 
of A*Prune remain exponential. Compared with some existing 
approximate algorithms [12] [14], the uncertainty of time and 
space complexities is the main drawback of A*Prune algorithm. 
In this section, we combine the A*Prune with any known e- 
approximate algorithms to give a Bounded A*Prune algorithm, 
BA*Prune, which can give either exact or c-approximate solu- 
tions to KMCSP problem in polynomial-time. However, the 
BA*Prune can only be applied to iZ < 2 cases due to the ap- 
plicability of known e-approximation algorithms. 

Let €-KMCSP be a known e-approximate algorithm to KM- 
CSP problem, BA*Prune will first run A*Prune for a bounded 
time T and bounded memory S (say the length of the candidate 
list Q), for example we can select S = and T = N^fxsec, 
If A*Prune is successful to give the exact solution within the 
bounded time and memory, we are done. Otherwise, BA*Prune 
executes the e-KMCSP to find an e-approximate solution. Since 
we can select S and T as polynomials of JV, then the BA*Prune 
is polynomial efficient algorithm if the e-KMCSP is polynomial 
algorithm. However the trade off for getting this efficiency is 
that BA*Prune is partially optimal rather than optimal: it is op- 
timal if it is run in the bounded time, otherwise is e-approximate. 
The completeness of BA*Prune is guaranteed by the complete- 
ness of e-KMCSP. All the other advantages of A*Prune are re- 
mained in BA*Prune, such as find K CSP paths and multiple 
constraints, as long as BA*Prune can be run in the bounded time 
and memory. 

VII. Comparison of Algorithms and Experimental 

Results 

In this section we compare the following algorithms in terms 
of CPU times: 

• A*Dijkstra: The A*Prune algorithm with Dijkstra distance, 
i.e., the A*Dijkstra algorithm presented in this paper. It gives K 
exact CSPs from one source node s to one target node t 
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TABLE I 

Algorithm running time comparison (p = 0.4. time unit: 

MICROSECOND) 



Input Data 


1-1 target 


1-N targets 


N 


D 


K 


A*Dijks 


L-Scale 


A*Unifo 


T-Scale 


100 


249 


1 

10 
100 


5357 
32773 
336367 


8443 


160734 


71020 


200 


241 


1 

10 
100 


20714 
54749 
365116 


#18420 


1084268 


1286781 


300 


238 


1 

10 
100 


38160 
186266 
1673267 


11705 


2994072 


1619317 


400 


236 


10 
100 


114843 
372441 
3291033 


# 85745 


4940668 


1197150 


500 


234 


1 

10 
100 


98889 
205030 
3670254 


#111651 


9548881 


2880185 



• L-Scale: The length scaling algorithm presented by R, Hassin 
[14]. It gives one CSP from one source node s to one target node 
t with length in the range of c-approximation. 

• A*Unifonn: The A*Prune algorithm with uniform distance, 
i.e., the A*Uniform algorithm presented in this paper. It gives 
K exact CSPs from one source node s to each of the N target 
nodes. 

• T-Scale: The delay scaling algorithm presented in [12]. It 
gives one CSP from one source node s to each of the N target 
nodes with delay in the range of e-approximation. 

We measured their execution times on a 450MHz SUN So- 
laris machine. For simplicity, we only tested for the case of two 
constraints, wq = wi is length and the W2 is delay. The error 
tolerance parameter e was set to 5% for the two e-approximate 
algorithms, L-Scale and T-Scale. 

For our experiments we used random graphs of a given node 
size N and connectivity p which is defined as the probability of 
a link existence between node pair (i, j). The link length, link 
delay and delay-limit D are also random numbers uniformly 
distributed in the open interval (0, 1000). The running times 
of the four algorithms for some randomly generated graphs are 
listed in Table I. To see how the running time varies with K, 
the A*Dijkstra algorithm is tested for K =1, 10 and 100 cases, 
while the other three algorithms are only tested for if = 1 case. 
The running time prefixed with a # indicates the non-optimal 
solution. 

These experimental results show that the A*Dijkstra is com- 
parable to L-Scale and the A*Uniform is comparable to T-Scale 
in running time. Considering A*Dijkstra and A*Uniform can 
find K exact CSP solutions, and can be applied to multiple con- 
straints and all the metrics can be float numbers, while the L- 
Scale and T-Scale can only give one e-approximate CSP solu- 
tion and can only be applied to two constraints with integer met- 
rics, we can say that our A*Dijkstra and A*Uniform have better 
performance on the average than the best known €-approximate 
algorithms, L-Scale and T-Scale. 



VIII. Conclusions 

A*Prune, an algorithm for finding K shortest paths subject to 
multiple constraints, has been presented. The algorithm grows 
a candidate path list, which contains the paths starting from the 
source node s and is initialized to the trivial path of (5, s). The 
candidate path list is ordered such that the path that most likely 
to project to a shortest feasible path is extended first. All the 
extended paths that are guaranteed to violate the constraints are 
pruned from the candidate list. Dijkstra's shortest path length 
can be used as lookahead feature in both the candidate path 
pruning and ordering processes. Experimental results show that 
A*Prune is comparable to the current best known polynomial- 
time e-approximate algorithms. BA^Prune, an algorithm com- 
bines the A*Prune with any other known c-approximate algo- 
rithm to give either optimal or e-approximate solution to the 
KMCSP problem in polynomial time, is also presented. 

Appendix 

We give proofs to the lemmas given in this paper. 
Lemma 1. Proof: Since p is a path from node s to node 
t, then 

T>{t,t,C(tyt)) = 0 
H(p) = W(p) 

Hence (12) is true. ■ 
Lemma 2. Proof: Suppose the edges on p is listed in 
the increasing order of number of hops from the start node s as 
{ei, 62, - . . , e/i}, where h is the number of edges on path p, then 

p = j>(s, s)eie2...efc 

Since path p{Sj s)ei is one step expanded path from path p{s, s), 
p{s, s)eie2 is two steps expanded path from path p{s, 5), etc., 
the path p is h steps expanded path from path p(s, s) . ■ 

Lemma 3. Proof: Assume path q is an inadmissible 
head path, i.e., 3r G (1,2,. ,R),Hr{q) > Cr- Let a: be 
an path expanded from q, i.e., 3y e P(V, V),a: = qy. Then 
3r € (1, 2, . . . , R), such that Hr{x) = Wr{qy) + Ar{d{y)) = 
Wr{q)HWriyHAr{d{y))) > Wriq)+Arid(q))) > Cr. This 
proves that path x must be an inadmissible head path. Lemma 3 
is true considering that x is any extended path of the inadmissi- 
ble head path q. ■ 

Lemma 4. Proof: As we said that because A*Prune 
expands paths in the admissible head path set P(s, V, H(p), C) 
in the increasing order of Ho{p), and all feasible paths also be- 
long to the admissible head path set, it must eventually expand to 
reach a feasible path. This is true unless there are infinitely many 
admissible head paths q with Ho{q 6 P(5, V^,H(p), C)) < 
Ho{p ^ P(s,i,H(p),C)). In the simple path search case, the 
admissible head path set must be limited, since its super set, the 
set of all simple paths in a graph is limited. However, in the com- 
plex path search case, loops with infinitely small length may re- 
sult in A*Prune searching infinitely many admissible head paths 
before finding a feasible solution. ■ 

Lemma 5. Proof: Given two paths u and v, both 
are feasible solutions of the KMCSP problem, and Wo{u) < 
Wo{v), Here Wo(p) is defined as the length of path p. A search 
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strategy for KMCSP problems is optimal if it finds u earlier than 

V. 

Assume A*Prune finds v earlier than u. The optimality of 
A*Prune can be proved by proving this assumption can not be 
true. We first decompose u into three parts: 

u = aeb 

where a 6 V,H(p),C), b e P(d(e),F), e e E. Then 
at the moment that A*Prune finds u, there exist paths a and b 
and an edge e, such that a is expanded but ae is not. Since u 
is a feasible path, its head path ae must be an admissible head 
path. Then ae is put into the candidate path list instead of being 
pruned when its parent path a is expanded. Since v is expanded 
but ae is not, then we must have: 

Ho{v)<Ho{ae) 



Hence: 



Wo{v) < Hoiae) 

< Woiae) Do{d{ae),t,C) 

< Wo{ae)-\-Woib) 



(16) 



This contradicts to the given condition. Then our assumption, 
A*Prune finds v earlier than u, can not be true. This proves the 
optimality of A*Prune. ■ 
Lfeiiima6. Proof: It is obvious that 



P(i,J,W(p),C(i,i))cP(i,j) 



then 



Vr € (0,1,.. .,/e) 

or equivalently, D(i,j,oo) < D(i, j, C(i, j)) must be re- 
mained true. This proves that the Dijkstra distance is admissible. 
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